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TECHNICAL NOTE 2726 


ON THE APPLICATION OF TRANSONIC SIMILARITY RULES 
By John R. Spreiter 


SUMMARY 


The transonic aerodynamic characteristics of wings of finite span 
axe discussed, from the point of view of a unified small perturbation 
theory for subsonic, transonic, and supersonic flows about thin wings. 
This approach avoids certain ambiguities which appear if one studies 
transonic flows by means of eq.iiations derived under the more restrictive 
assumption that the local velocities are everywhere close to sonic veloc- 
ity. The relation between the two methods of analysis of transonic flow 
is examined, the similarity rules and known solutions of transonic flow 
theory axe reviewed, and the asyn 5 >totic behavior of the lift, drag, and 
pitching-moment characteristics of wings of large and small aspect ratio 
is discussed. It is shown that certain methods of data presentation are 
advantageous for the effective display of these characteristics. 


INTRODUCTION 


Th.e' small perturbation potential theory of transonic flow proposed 
apparently independently by Oswatitsch and Wieghaxdt, Busemann 
and Guderley, von Karman (references 1 through 6), and others is now 
supplying a rapidly increasing fund of information regarding transonic 
flow about aerodynamic shapes. Solutions have been given recently for 
the flow around symmetrical nonlifting airfoils at both subsonic and 
supersonic speeds in papers by Guderley and Yoshihaxa, Vincent! 
and Wagoner, Cole, Trilling, Oswatitsch, GiiU-strand (references 7 
through 14), and others. In the application of these- results to specific 
examples, two items of theoretical interest have been noted (see, in 
particular, references 8, 15, and 16); (a) The theoretical results 

appear to be applicable at Mach n umb ers fax removed from 1 even though, 
in most cases, the results have been obtained from equations valid only 
in the immediate neighborhood of sonic speed, (b) In the application 
of theoretical results to specific exan 5 >les at Mach numbers other thart 1, 
it has been noted that certain ambiguities exist in the theoretical 
determination of aerodynamic quantities . It is one of the purposes of 
this report to investigate these two points. This is accomplished by 
examining transonic flow from the point of view of equations 



2 


MCA TW 2726 


that are valid throughout the Mach number range rather than only in the 
neighborhood of sonic speed. Such an approach emphasizes the relation 
between the roles of linear theory and of nonlinear theory in the tran- 
sonic range. 

The similarity rules provided by the theory (references 5 ^ 6 , 
and 17 through 20 ) have also proved to be useful in the correlation and 
interpretation of experimental data. It is with the latter aspect of 
the transonic-flow problem that the present paper is primarily concerned. 
In this paper, the similarity rules and their application to the specific 
problem of concise presentation of lift, drag, and pitching-moment char- 
acteristics of wings are given in detail. The known solutions of two- 
dimensional transonic flow are reviewed and the asymptotic behavior of 
the aerodynamic characteristics of wings of large and small aspect ratios 
is examined. It is shown that certain methods of 'data presentation are 
advantageous for displaying these characteristics. 


SYMBOLS 


A aspect ratio 

A [( 7 +l)(t/c)]^/^A 

a speed of sound 

Sq speed of soxind in the free stream 

a* critical, speed of sound 

b wing semispan 

Cp drag coefficient 

Cp^ drag coefficient of symmetrical nonlifting wings 

ACp contribution to drgig coefficient due to lift 

(ACp/a^) [( 7 +l)(t/c) ]^/®[ACp/a^] 

Cl lift coefficient 

(C^) f(7+l)(t/c)]''/®[CL/d] 

Cm pitching -moment coefficient 

( Cm/a) [ ( 7 + 1 ) ( t/c) ] [Cm/a ] 


pressure coefficient 
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C.P. 

c 

°<io 

Aca 

(Acd/a^) 

Cl 

(cj/a) 


center-of -pressure function 
wing chord 

section drag coefficient of symmetrical nonlifting airfoils 

[(r+i)^/7(t/c)=/=l 

contribution to section drag coefficient due to lift 
[ ( 7+1) ( t/c) ] ^/°[Acd/a^] 
section lift coefficient 
[(7+l)(t/c)]'-/^[cz/a] 


D 

Do 

Da 

do 

dA 

L 

1 

M 

m 

Mo 

• 

P 

s 

t 

Uo 

x,y,z 


^c.p. 

(z/c) 


drag function 

draLg function for symmetrical nonlifting wings 
drag due to lift function 

section drag function for symmetrical nonlifting airfoils 

section drag due to lift function 

lift function 

section lift function 

pitching -moment function. 

section pitching-moment function 

free-stream Mach number 

pressure function 

stretching factors defined in equation (B 7 ) 
maximum thickness of wing 
free-stream velocity 

I 

Cartesian coordinates where x extends in the direction of 
the free-stream velocity 

distance from wing leading edge to center of pressure 
ordinates of wing profiles in fractions of chord 
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a angle of attack 

a a/(t/c) 

r 7 + 1 

7 ratio of specific heats, for air 7 = 1.4 

X arbitrary constant 

lo (Mo^-l)/[(7+l)(t/c)]^/3 ^ 

T ordinate-amplitude parameter 

4 velocity potential 

<p perturbation velocity potential 

Subscripts 

I values given by linear theory 

¥ conditions at the wing surface 

FUWDAMEiraAL CONCEPTS 
Basic Equations 


The quasi-linear partial, differential equation satisfied by the 
velocity potential. $ of steady isentroplc flow of a perfect Inviscid 
gas can be expressed in the form 



■'xx 


4 2 

■ y 

■ a^ 

J ®yy + 


'tyi'z 

a2 

0 

1 

ro 

-4^ Ozx = 0 

a2 


zz 


_ 2 


®x®y 


4 >xy - 


( 1 ) 


where the subscript notation is used to indlca;te differentiation 
and a is the local speed of sound given by the relation 



a2 = a 2 _ Zzi 

o 2 


( 2 ) 
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In "tills la"t"ter equa"tion Uq and Sq are, respectively, "the veloci"ty and 
speed of sound in the free stream and 7 is the ratio of specific heats 
(for air 7=1 A). 

Introducing "the perturbation velocity potential 9, "where 

cp = -DqX + 4> ( 3 ) 


it is possible to eiqpress equation (l) in terms of the deri"vatives of 9 
as follows; 


*Pyy 


(7"H)Uo<I^ 9x^ + (iy^+ 9z^)] + 


1-Mo^J q>xifKPyy+<Pzz=< ^ < 


(-7-1)Uq9x+ ') ^ 


J" 


w 


2 ^ ,1’. [«o + fx] + 





9y 9z 


J 


If it is assumed that al 1 perturbation velocities and perturbation 
velocity gradients (represented by first and second derivatives, respec- 
tively, of 9 ) are httih. 1 1 and that only "the first-order terms in s mal l 
quantities need be retained, equation (4) simplifies to the well-kno"wn 
Prandtl-Glauert equation of linear theory 


( 1-Mn^ ^9 +9 + 

I ^ J ^xx yy zz 



(5) 


\ 
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Z 



(a) 


where the free-stream velocity is 
directed along the positive x axis 
as shown in sketch (a) and where Mq 
is the Mach number of the free stream. 

It is well-known that equation (5) 
leads to iiseful results in the study of 
subsonic and si^ersonic flows about 
thin wings and slender bodies but that 
it is incapable, in general, of treat- 
ing transonic flows. The failure of 
linear theory in the transonic range is 
evidenced by the calculated value of cpx 
growing to such magnitude that it can 
no longer be regarded as a small quan- 
tity when conpared with Uq. 


Second-order theory for thin wings would involve solution of the 
equation 


^1-Mq^^ 9xx T’x ?xx + ^ 


( 9y q>xy + *Pxz) 
Uq 


(6) 


Actually, we are interested in retaining higher-order terms only to 
the extent that is necessary to allow the study of transonic flow. 
Examination of the known characteristics of transonic flow fields indi- 
cates that the first term on the right can often become of importance 
and should be retained. The remainder of the terms on the right can 
never become large for treinsonic flows about thin wings at small angles 
of attack and can be safely disregarded. Furthermore, since the right- 
hand side is merely an approximation to allow the treatment of transonic 
flows and rapidly diminishes in magnitude as Mq departs from unity, 
the eq^uatlon may be further simplified without much, if any, loss in 
accvuracy by setting Mq = 1 in the coefficient of the term on the right, 
(if one does not wish to make this added approximation, the results of 
this paper should be altered by replacing 7 + 1 with Mo^(7+l) 



MCA TN 2726 


7 


wherever it occurs.) The slr^jlified equation^ is 

(l-Mo^) «Pxx + "Pyy + «Pzz = <Px 'Pxx (7) 

In addition to satisfying the differential equation, the perturba- 
tion potential must provide flows conipatihle with the following physical 
requirements: (a) The flow must he uniform far aheeid of the wing and 

(h) the flow must he tangential to the wing surface. Therefore, the 
following boundary conditions are to he specified for the perturbation 
potential; 

at X = - 00 


('I’xio = ' Wo - ° 


(8) 


at the wing surface, W 


_1 

Uo 


(cPjj) = (^Z/Sx) 
W 


( 9 )- 


where (Sz/^x) refers to the local slope in the x direction of the wing 
surface. Furthermore, it is consistent with the assumption of small 
disturbances to satisfy the second boundary condition on the two sides 
of the xy plane rather than on the wing surface. Equation (9) is 
therefore replaced by 


i ^ ^(f. f) 


(10) 


^Although equation (7) is valid throu^out the Mach number range, it is 
not the appropriate equation for the treatment of the "pseudo tran- 
sonic" flow fields which are to be found aroimd three-dimensional 
swept wings (consider, for instance, an infinite yawed wing) throu^ 
a limited range of supersonic Mach numbers. Since these flows are 
characterized by shock waves, standing essentially normal to 
the xy plane but oblique to the free-stream direction, equation ( 6 ) 
can be simplified only to the following: 

('l- v) ' ‘Pxx 'Pyy + f zz = 'Px 1>xx + 

T “Px ‘Pyy + ^ fy fxy 

The discussion of these problems is outside the scope of the present 
report . 
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where the shape of the ^^Ing profile is given hy 

(z/c) = Tf(x/c, y/b) ( 11 ) 

where f(x/c, y/b) represents the ordinate-distribution function and t 
is an ordinate -amplitude parameter.^ In order to obtain uniq_ue and physi- 
cally in^jortant solutions, it is necessary t& assume the Kutta condition 
that the flow leaves all subsonic trailing edges smoothly. To be com- 
plete, consideration should also be given to the conditions of transition 
through shock waves. This point has been discussed by Guderley (refer- 
ence 4 ) and hy Vincenti and Wagoner (reference 9) • They indicated that 
(a) eq^uation (7) is valid even if shock waves are present, and (b) the 
shock conditions aigree with the similarity rules. In' this report, the 
similarity rules are derived in appendix B and their conq)atibility with 
the shock relations is demonstrated in appendix C. 

Upon solving the above boundary-value problem for the potential, one 
may determine the pressure coefficient by means of the formula 



( 12 ) 


Althoiagh Osira-titsch, Berndt, and GuUstrand (references 12 , 13 ; 1 ^; 
19, and 20) have previously investigated transonic-flow phenomenon by 
means of equations derived by assuming that all velocities are small per- 
turbations around the free-stream velocity Uq as is done above, most 
other workers have used equations derived under the more restrictive 
assunq)tions that all velocities are small perturbations around the criti- 
cal velocity of sound a*. In the latter scheme, the perturbation poten- 
tial is defined by (see, for instance, reference 6 or I8) 


cp' = -a*x + 3 > (13) 

and the resulting differential equation for cp' is 




+ 'Pzz = ^ 'Pv 9; 
• a* 


XX 


The corresponding boundary conditions are specified as follows: 




^Wote that, in general, a variation of t represents a simultaneous 
change of the thickness ratio, camber, and angle of attack. In the 
special case of nonlifting wing .having symmetrical sections, t is 
proportional to the. thickness ratio; for inclined flat-plate wings of 
vanishing thickness, t is proportional to the angle of attack. 
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at X = - 00 


(<f^x >o = Uo-a* fw _ _a_ 




= 0 


(15) 


at the wing surface . 


(<P’J 


z'z=o 


a^T 


^(x/c) 




(16) 


where the shape of the wing profile is still given by eguatlon (U) . 
The equation for the pressure coefficients is approximated similarly, 
thus. 


(IT) 

The two statements of the equations for transonic flow are clearly 
identical at a Mach number of unity. Although the derivation of the 
a* equations requires that the free-stream Mach number be very close to 
unity, these equations have been used with good success by a number of 
authors to calculate the aerodynamic forces on airfoils at Mach numbers 
considerably removed from unity. In so doing, it has been suggested that 
it might be preferable to vise more accurate relations for the pressure 
coefficient or the boundary conditions; for instance, it has been suggested 
that a* be replaced with Uq in the equation for Cp’. This matter 
has been discussed at length in references 8 , I5, aui I6. Since no 
restriction requiring the Mach number to be near unity is made in the 
Uq analysis, it is informative to examine the relation between the 
results of the a* and the Uq analyses. This is done in appendix A. 

It is found that the a* analysis, if performed in a cong)letely con- 
sistent manner using equations (13) throu^ (l 7 )# yields values for Cp 
that are identical to those given by the more general Uq analysis. 

This somewhat paradoxical result is achieved through the action of a 
number of compensating effects and only applies to the pressures and the 
forces and moments derivable therefrom. It should be noted, in particu- 
lar, that the values of the local velocities and Mach numbers provided 
by the a* analysis for flows having free-stream Mach numbers other than 
unity are in error. Throughout the remainder of this report, the di’s- 
cussion will be based on the Uq analysis. 

It is ln5)ortant to recognize that wing theory based on equation (7) 
is valid for all Mach numbers below the hypersonic range. At subsonic 
and supersonic speeds, equation (7) is of the same order of accuracy as 
the Prandtl-Glauert equation of linear theory (equation (5)) although more 
difficult to solve. At Mg = 1 , equation (7) is identical with equa- 
tion (1^), now widely used in the study of transonic-flow problems. 


a — 


q>’ - (cp;) 

X X Q 
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On the other hand, there is no a priori method for determining 
•whether or not a solu'tion of the equations of linear theoiy ■will he 
'valid in the transonic range. One can only decide hy solving the prob- 
lem under the assumptions of linear theory eind' then inspecting the magni- 
tudes of the terms, particularly of qpjj, to see whether or not they can 
be regarded as small quantities. If the terms are sufficiently small, 
the linear -theory solution is presumed valid even though the Mach number 
may be near unity. Linearized-theory solutions have been obtained for a 
great number of practical •wing problems and their behavior in the tran- 
sonic reinge is now well known. To re'view briefly: For unswept wings of 

infinite span, linear theory Indicates that the ma gnitude of cp^ on the 
surface of a given airfoil is proportional to l// |l - Mq 2| ; consequently 
cPx approaches infinity as Mq approaches unity and the theory is 
clearly Inapplicable. For •wings of finite span, however, the perturba- 
tion velocities may be large or small at sonic velocity, depending on 
the particular problem as discussed in detail in reference 21 . Specifi- 
cally, for three-dimensional lifting surfaces of zero thickness, the 
velocities remain finite everywhere except at the leading edges, their 
magnitudes generally Increasing ■with Increasing aspect ratio and angle 
of attack. For -wings of nonzero thickness, however, cp^j. generally 
becomes large logarithmically as 1 - Mq^ approaches zero; consequently, 
linear theory is inapplicable -within some Mach number range s'urrounding 
unity. 

Summarizing, linear theory is applicable to lifting surfaces of 
small or moderate aspect ratio at all transonic speeds, but fails for 
wings of finite thickness wi-thin a range of Mach number sxorrounding 
unity. The range of inapplicability diminishes to zero as the aspect 
ratio, thickness ratio, end angle of attack of the wing tend to zero. 

In treating transonic flows for which linear theory is applicable, 
it is often advantageous to consider the special case of sonic flow 
(Mq = 1) separately. Equation (5) for the perturbation potential then 
reduces to a particularly sinqile form 


'Pyy 


q> 


zz 


= 0 


(18) 


Solutions of equation (I8), in conjunction with the boundary conditions 
given by equations ( 8 ) and (lO), are identical to those of linear theory 
foxmd by solving equation (5) and subsequently setting Mq = 1, but can 
be ob-tained with much less effort. Since, in addition, the results of 
this sin5)le theory, now generally knovm as slender -wing theory, are also 
applicable to low-aspect -ratio lifting surfaces throughout the 
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entire Mach number range,® a considerable number of solutions of slender- 
wlng theory have been presented in the._last few years. (See refer- 
ences 22 , 23 > and 2 ^ 4 .) These results are, of course, applicable to flows 
at Mq = 1 to exactly the same extent as the results of linear theory. 


Similiarity Rules 


In reference 6 , von Karman derived similarity rules for the pressure 
distribution, lift, drag, and pitching moment of airfoils in treinsonic 
flow using equations (l 3 ) through ( 17 )* The same equations were used in 
reference I8 to determine the transonic similarity rules for wings of 
finite span. The corresponding similarity rules of linearized subsonic 
and supersonic wing theory were also derived and compared with the tran- 
sonic similarity rules in the latter reference. It was shown that the 
similarity rules of linear theory contain an arbitrary parameter and can 
be expressed in many forms, one of which coincides with the similarity 
rules of transonic flow. 


It follows from appendix A and is demonstrated in detail in appen- 
dix B that the similarity rules derived from equation (7) are Identical 
to those previously given in references 6 and I8. The new derivation, 
however, possesses the advantage , of being based on a single statement of 
the problem of wing theory that is uniformly valid at subsonic, tran- 
sonic, and supersonic speeds. The similarity rules for Cp, Cj^, Cjjj, 
and Cp were stated in reference 18 to be 


Cp 


tS/® ^/l-Mog 

(7+1)^/® ^ [(7+1) t ]^/®^ 


yi - Mo^ A; 


2 E y 

c' b 


(19) 


Cl 


t 2/3 

( 777 ^ 


yi-Mpg 

[( 7 + 1 )t]^/®' 


yi - Mo^ A 


( 20 ) 


®Thls dual role of slender-wing theory stems from the two mys that one 
can reduce equation (5) to equation (18). One can neglect the term 
(1 - Mo^) cPxx conqjarison with q)yy and cp^z either because cpxx 
is small, as may be the case with law-aspect-ratio wings or slender 
bodies, or because (l - Mq^) 1b zero, provided cpxx does not become 
very large in comparison with the other velocity gradients T>yy 
and cpjjz* application of this theory to slender wings 

antecedes the application to flows with sonic free-stream velocity, 
hence the name slender-wing theory. 
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t2/3 

' (r+i)V3 “ 


t5/3 

Cd = D 


i/l-Mo" 


[(r+i)T]V3 




1 - Mo A 






, yi - Mo^ A 


where the geometry of related wings is given by equation (u): 

(z/c) = f f(x/c, y/b) 


( 21 ) 


( 22 ) 


Equations (19) through (22) are functional equations. For exaniple, 
equation (19) is to be interpreted as stating that the pressure coef- 
ficient Cp is equal to 7 +I) times some fimction P of a 

number of specified parameters. The foregoing equations have been 
written for fioi'TS where Mq< 1. If Mq^I, the radical J \ - Mq^ 

should be replaced with a/Mq^ - 1. The functions P, L, M, and D eire 
different, however, for subsonic and supersonic flow. Consequently, 
subsonic flows may be related to other subsonic flows by the similarity 
rules, but not to supersonic flows, and conversely. 


It is inportqnt to recognize that the similarity parameters may be 
conibined or regroT:ped in ainy manner whatsoever, provided the same number 
of independent parameters is always retained. For instance, in much of 
w hat follow s, it will be found desirable to use t he square of 
4 1 - Mo^/[(7+1)t]^/® and to replace a/i - A with a new parame- 
ter [(7+1 )t]^/®A obtained by dividing 4 1 - Mq^ A by 
a/ 1 - Mo^/[( 7+1) In terms of these parameters, the similarity 

miles are 


pi 

Mo^-l 

^ (7+1)V3 ^ 

[(7+1)t]2/3 

t2/3 

*^7 - L 

Mo^-l 


[(7+1)t]2/3 

c = M 

Mo2-l 

(7-M)"/" , ' 

[(7+1)t]2/3 


[(7+1)t]V=’A; 2, I 
c b 

[(7+1)t]^/3 ^ 
[(7+1)t]V3 a] 


( 23 ) 

(24) 

( 25 ) 
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Cd = 


t5/3 


(y+l)l/3. 


D 


MqS-I 


[(7+1)t]2/3 


[(7+1)t]V3 a 


(26) 


wherein the geometry of related wings is again given by eq.uation (u) . 


Slope of Presaxire Curve at Mq = 1 


Liepmann and Bryson (reference I 5 ) have recently detennined the 
slope of the Cp versus Mq cnrve at Mq = 1 by means of the foUo^dJig 
simple and intiiitive considerations. It is a we3JL-known fact that, at 
slightly supersonic Mach nuuibers, the detached bow wave is far away 
from the airfoil and ne^ly normal. It is also well known that the Mach 
number downstream of a weeik no rmal, shock is as much below unity as the 
Mach number XQ)stream is above unity. Consequently, the pressure or Mach 
number distribution on the airfoil should be independent of Mach nuuiber 
in the neighborhood of Mq = 1. The slope of the Cp versus Mq curve 
at Mq = 1 can then be found by simple and direct means and is the fol- 
lowing for thin airfoils: 


( = ... • 


(27) 


Vincent! and Wagoner (reference 8) have found by means of similar 
considerations that a more exact relation is given by 


^ 2_ 

W/mo=x ^ ^ 

Since the above derivation is based to a certain extent on physical 
reasoning which may be either more or less exact than small perturbation 
transonic theory, it is of interest from the present point of view to 
determine the equivalent resxilt from the model of transonic flow provided 
by equation (7). The result, obtained by a similarity type of analysis 
and presented in detail in appendix D, is just that which would be found 
if one actually solved equation (7). . 

It is found that a solution for Mach numbers slightly less than 
unity satisfies the differential, equation and the boundary conditions at 
the wing surface for flows with Mach numbers slightly greater than unity. 
However, the velocity perturbations do not go to zero Infinitely fax 
Eihead of the wing, but gg instead to the value corresponding to that 
associated with a nonnal shock ’t^ave. If it can be assumed that the bow 
wave approaches a normal shock wave standing infinitely far ahead of the 
wing at a sufficiently rapid rate as the free-stream Mach number 
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approaches unity, the slope of the Cp versus Mq curve at Mq = 1 is as 
given hy equation ( 27) • Whether or not this is a permissible assun^ition 
for any given case still remains an unanswered question. Intuitive con- 
siderations suggest that the results are probably applicable to symmet- 
rical airfoils at zero or Infinitesimal angles of attack but not to air- 
foils at larger angles of attack or to wings of finite span. 


APPLICATIONS 

Fundamental Hypotheses and Principles 


The remainder of this report is princip ally concerned with the 
deduction of the qualitative, and to some extent ^quantitative, character- 
istics of thin wings in transonic flow by means of simple logical con- 
siderations based primarily on the similarity rules together with the 
following hypotheses: 

(a) Nonlinear theory based on equation (j) is applicable to all 

problems . 

(b) Linear theory based on eqixation ( 5 ) is valid for all wings at 

Mach numbers either appreciably below or above unity. 

(c) Linear theory is valid at all, Madh numbers, except possibly 

very near unity, for wings of small aspect ratio. 

(d) The differential pressxires between the upper and lower surfaces 

of a wing having symmetricaJL airfoil sections are propor- 
tional to the angle of attack for at least a small range of 
angles about zero. 

(e) The slope of Cp versus Mq at Mq = 1, defined by equa- 

tion (2j), is applicable at least to symmetrical airfoil 
sections at zero or infinitesimal angles of attack. 

The consequences of the foregoing statements will be consistently 
pursued in the following sections in the discussion of the aerodynamic 
characteristics of airfoils and complete wings. Throughout, the analysis 
will be restricted to wings having symmetrical profiles. The decision to 
ignore the influence of camber is based not only on the desire for sim- 
plicity but on the fact that symmetrical airfoils appear experimentally 
to have superior aerodynamic characteristics in the transonic range. 
Whenever specific results are to be used to illustrate the statements, 
they will nearly always be for symmetrical-wedge or double-wedge profiles 
pTifl for wings of triangular plan form. This choice is dictated by the 
present availability of theoretical results. 
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The basic principle in the following analysis is to express the 
similarity rules in such forms that the lift, pitching -moment, and drag 
coefficients can be studied for limiting values of the parameters with 
no chance for ambiguity due to indeterminate forms. In this respect, 
the statement of the similarity rules provided by equations ( 23 ) 
through ( 26 ) will be found particulaxly useful. 

The similarity nles thus formulated are totally equivalent to those 
given by equations ( 19 ) through (22) but possess three outstanding advan- 
tages: 

(a) The indeterminacy at Mq = 1 resulting from two parameters 

simultaneously vanishing is removed. 

(b) The squaring of the first parameter avoids the necessity of 

changing parameters as sonic speed is passed. 

• « 

(c) The use of t he parameter [( 7 + 1 )t]^/® A rather than 

a aids in distinguishing the regimes in which 
linear theory is applicable in the transonic range from 
those in which nonlinear theory must be used. Thus as 
[(7+1 )t]1/3 a approaches zero, linear theory is always 
applicable provided, in some cases, that Mq is not pre- 
cisely equal to unity. On the other hand, as E(7+1 )t]^/® A 
becomes large, linear theory is not applicable in the tran- 
sonic range and nonlinear theory must be used. 


Pressure Drag of Symmetrical Nonlifting Wings 


The similarity rule for the pressure drag coefficient of symmetrical 
nonlifting wings having profiles given by 


(z/c) = (t/c) f(x/c, y/b) 


(29) 


is obtained from equation (26) by Identifying T with t/c, the thickness 


ratio. 


»0 (^+ 1 ) 1/3 


D. 


Mo^-1 


[(7+l)(t/c)]^/=^ 

r\ (e X) 


> [ (7+1) (t/c) A 


( 30 ) 
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Therefore, dreig resulhs for Bynnnfitrical nonlifting wings should he 

presented hy glotting the variation with and A of a genersilized drag 

coefficient Cn defined as follows: 

•^o 


CDo 


^ (y+l )^/3 
“ (t/c )=/3 


CDo = Do do. 


A) 


(31) 


At Mach numbers stifficiently removed from unity for linear theory 
to apply, Cpo must he independent of 7 since 7 does not appear in 
either the differential equation or boundary conditions of linear theory. 
Therefore, 



(t/c)=/=‘ 

lv-ll 

- 1 /a 

D 



( 7 + 1 )^/® 

[( 7 + 1 ) (t/c) 

h>oj 


[(7+l)(t/c)]2/3 


[(7+l)(t/c)]^/® A 


■ (Vc)^ Dq 

y (mo^-iI ^ 


A^ 


(32) 


where it should he recalled that Dqj is a different function for suh- 
sonlc ajad supersonic flow. Equation ( 32 ) is equivalent to the extended 
Prandtl-Glauert rule. For suhsonic flow, D'Almbert’s paradox requires 
that the drag he zero; therefore, for eill wings, 

(Cd„)j=0, (33) 

Mo<i go<o 

For supersonic flow, wave draig exists idilch depends on / Mo^ - 1 A 
as well as on the plan form and airfoil section. The general functional 
relation for the drag coefficient of a family of affinely related wings 
at zero ang le of attack, as given hy linear theory, is 

(OBo), = (7^ a), (|7''"a) (34) 

Mo>i V Mo -1 \ / 

Wings of infinite aspect ratio .- For wings of infinite span (or 
airfoils), equation (32), representing the functional relation of linear 
theory for the drag coefficient, reduces to the following: 
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(cdo), =-===>< const., (35) 

Mo>i lo>o ^ 

where the value of the constant depends on the shape of the airfoil. 
Numerous experimental, data show variations consistent with equation ^ 35 ) 
at Mach numbers greater than about that of shock attachment. At Mach 
nimbers closer to unity, however, the theoretical values provided by 
this eq.uation are unreliable. It is evident from inspection of the 
results of linear theory itself that such a failure occurs, since the 
perturbation velocities assumed to be small in the derivation of the 
eq.uations are found to become infinitely large as the Mach number 
approaches unity. It is apparent, therefore, that it is necessary to 
resort to nonlinear theory for the calculation of the drag of airfoils 
in the transonic speed range. 

A similarity inle for the section drag coefficient of symmetrical 
nonllfting airfoils which is veilid throughout the Mach number range may 
be obtained from equation ( 3 I) by setting A equal to infinity. 

^ “ Do (So. ”) = do (So) (36) 

At a Mach number of unity, the similarity parameter vanishes 

and the function do (Iq) Is si constant. 

(^) = do (0) = const., (cdo) = X const. (37) 

go=o Mo=i (t'+I) 

indicating that the section drag coefficients of affinely related air- 
foils are proportional to the five-thirds power of their thickness 
ratios. If hypothesis (e) is accepted> the variation of c^^ with Mo 
at Mo = 1 is found to be zero for complete airfoils 



Since calculations have been made of the drag in transonic flow of 
sin5)le symmetrical sections at zero angle of attack, it is not necessary 
to speculate further regarding the variation of c^^ with ^o* 
present, however, the profile for which the most complete information 
exists is not a conrplete airfoil but is a single-wedge section followed 
by a straight section extending infinitely far downstream. In accord 
■(•rLth some of the original papers op this subject, the single-wedge section 
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is considered as the front half of a symmetrical douhle-wedge airfoil 
having a chord c. The value of t/c in c^ and Iq is defined accord- 
ingly. Solutions for this section obtained using the nonlinear smaUL 
perturbation theory have been given for flows having subsonic, sonic, and 
supersonic free-stream velocities, respectively, by Cole (reference 10 ), 
Guderley and Yoshihara (reference 7 ), and Vincenti and Wagoner (refer- 
ences 8 and 9 ) • The linear-theory solution for pure supersonic flows 
has been given by Ackeret (reference 25) . The slope of the drag curve 
at Mq = 1 is no longer zero as indicated for the complete airfoil by 
equation (38) but takes on a positive value given originally by 
Liepmann and Bryson (reference 15) and readily derivable from eq.ua- 
tions (27) and (38). 



All of these results are combined on a single graph in sketch (b) . It 
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may be seen that the preceding remarks concerning the relation of linear 
theory and nonlinear theory and the slope of the drag curve at a Mach 
number of unity are verified by this cong)arison. 

An indication of the accuracy of the theory is provided in 
sketch (c), which shows the theoretical curve of sketch (b) together 
with corresponding results obtained from wind-tunnel experiments by 
Liep mann and Bryson (references I5 and I6) . The vertical lines through 
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the experimental points indicate estimated accuracy of the data. Three 
single-wedge models having semiapex angles of ^-5°^ 10° were 

used in the tests. This figure illustrates the degree to which the 
s im ilarity rules axe able to reduce data from a family of profiles having 
different thickness ratios to essentially a single curve. 

The double-wedge airfoil has also been treated theoretically through- 
out the entire Mach number range. Solutions for Mq < 1 have been given' 
by Trilling (reference U), for Mq = 1 by Guderley and Yoshihara 
(reference 7)^ and for Mq > 1 by Vincent! and Wagoner (references 8 
and 9) • The results of their calculations are shown in sketch ( d) . 



2 
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For Mq < 1 , Oswatitsch has given approximate solutions for the 
pressure distribution on symmetrical biconvex airfoils (references 12 
and 13) and for the pressure distribution and dreig on MCA four-digit 
symmetrical airfoils (reference 13 ) • This work has recently been 
extended to several MCA 6 -series airfoils by Gullstrand (reference 1 ^<-) . 
Their drag results axe generaULy similar to those indicated at corre- 
sponding Mach numbers for the double-wedge section in sketch (d). 

Wings of finite aspect ratio .- The similarity rules for wings of 
finite aspect ratio are given by equations (31) and ( 32 ). Although no 
essential simplification of the rules occurs for wings of small aspect 
ratio, the range of applicability of linear theory increases as the 
aspect ratio decreases. This point can be illustrated by considering 
the results provided by linear theory in a specific case. A good example 
to select for this purpose is that of the drag in supersonic flow of a 
triangular Wing -irith symmetrical double-wedge airfoils. (See refer- 
ence 26.) This particiilar choice tos made for the foUoi'dng reasons: 

(a) Solutions are known for at 1 supersonic Mach numbers; (b) the double- 
wedge airfoil discussed in the preceding sections corresponds to the 
limiting case of the Tfing of very great aspect ratio. The dreig results 
provided for this ■tring by linear theory are presented in sketch (e). 




The results of sketch (e) are presented in a different manner in 
sketch (f) wherein Cp^ is plotted as a function of various 

values of A as suggested by equation ( 3 l). For purposes of comparison, 
the curves for the drag of airfoils (A=a>) computed by both linear and 
nonlinear theory are also included on the graph. As noted in the pre- 
ceding section on airfoils, comparison of the results of linear theory 
inLth those of nonlinear theory for wings of A = 00 shows that good 
agreement exists for larger ^ q , but that at smaller ^q, the values 
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predicted by linear theory become too large. For wings of finite A, 
only the results of linear theory are available. They also exhibit the 
trend of indicating infinite drag as |o approaches zero; however, the 
range of in which the value of Cp^ is excessive is much less than 

is the case for A = <». In general, 5p^ of wings of small aspect ratio 
diminishes with decreasing aspect ratio. 

The drag results of sketch (e) are presented in still another form 


in sketch (g) in which is plotted the 
various values of 5 q. The princi- 
pal merit of this method of plot- 
ting is that it aids in distinguish- 
ing the region -^diere nonlinear theory 
must be used from that where linear 
theory may be useful. ‘ Thus the txro- 
dimensional nonlinear theory results 
appear on the right of the graph 
corresponding to large 2 l, whereas 
the three-dimensional linear theory_ 
results appear- on the left for 
small A, The filling in of the 
remainder of the graph requires 
either the solution of the equations 
of three-dimensional nonlinear -i-ring 
theory or the use of experimental 
data. It should eigain be noted that 
the present drag considerations apply 
only to the pressure drag. Before 
plotting experimental results in the 
manner indicated, it is necessary to 
first subtract the friction drag. 


variation of CDq with A for 



Lift 


Equation (24) indicates that the similarity rule for the lift 
coefficient of wings having profiles given by 

(z/c) = Tf(x/c, y/b) 

is 


(^+ 1 ) 1/3 


[ (7+1)t]^/® 
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This eq.uation has been the soxirce of sonie confusion due to the multiple 
role that t plays in determining the thickness ratio, camber, and angle 
of attack. A more explicit statement of the similarity rule that is 
useful for a1 1 -wrings having symmetrical profiles of nonzero thickness is 
given by the foUo^vLng pair of equations: 


(z/c) = (t/c) f(x/c, y/b) 


(40) 


Cl = 


(7+1) 


L’ 


^ ^ ~ ^ , [( 7 +l)(t/c)]^/° A, 

[(7+l)(t/c)]2/3^ t/c 

(7+1) 


(41) 


where the primes serve as a reminder that L and L are different 
functions of the parameters indicated. If hypothesis (d) is accepted, 
Cl varies linearly with a for at least small singles of attack. It 
is advantageous , therefore, to consider the lift ratio Cl/cc rather 
than Cp Eilone, thereby minimizing the influence of a. 


Cl ^ 1 (t/c)^/^ 

O' ~ O' (7+1)^/° 


L'Uq, a, a) = 


ttt l"(6 . s:, S) 

[( 7 +l)(t/c)]^/° ° 


(42) 


Therefore, lift results may be presented by plotting__the variation 
with 5o» A, and a of a generalized lift ratio Cl/o. defined as follows 
( the primes on L have been omitted for simplicity) : 


ClT^ = [( 7 + 1 ) (t/c) (CL/a) = L(lo^ “) ( 43 ) 

Equation (l<-3) shows that Cl/<x depends upon three parameters, one 
more than the number which can readily be treated on a siB 5 >le plot. 
Simplification can be gained, of course, by holding one of the parameters 
constant for an entire graph. Results so presented are psirticuleirly 
interesting for |q= 0, (Mq=1); S=“, (A=»)j and a=0, (cx=0) . The latter 
scheme is especial 1 y good since experiments indicate that lift curves of 
wings are often relatively straight lines at a1 1 Mach numbers. The 
values of Cp/a at a = 0 might, therefore, be expected to be good indi- 
cations of the acttial values for other a. The appropriate similarity 
nale may then be written 

(cij) -L(S„, A, 0) =L„(£c^ it) 


(lA) 
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For cases where linear theory applies (hypotheses (b) and (c)), two 

statements can he made immediately which provide further information 

about CL/a: (a) {CjJoi) must be independent of y, and (b) (C£/a) 

b b 

must be Independent of a by virtue of the superposition principle of 
linear theory. Therefore, following the procedure used in equation (32) 
gives 






where again L ^ is a different function for subsonic and supersonic flow. 


Wings of infinite aspect ratio .- J’or wings of infinite aspect ratio, 
the functional relation of linear theory for the lift ratio given by 
equation (45) reduces to 



Solutions of the equations of linear theory show that the value of the 
constant is 2it for subsonic flow and four for supersonic flow. Exami- 
nation of these results indicates that they are valid at Mach numbers 
appreciably less than or greater than unity, but are invalid for Mach num- 
bers near 1. 

A similarity rule for the section lift coefficients of a family of 
affinely related symmetrical, airfoils which is valid throughout the Mach 
number range may be obtained from equation (43) by setting A = “>. 

(cj/a) = L(5 q, o) = 1 (Iq, a) (.4?) 


At a Mach mamber of unity, |q is identically zero and the 
expression for the lift ratio becomes 


(c^/a) =1 ( 0 , a), 
^0=0 


(cz/a.) 

Mo=l 


1 

[(7+l)(t/c)]V3 





Equation (48), when considered together with hypothesis (d), indicates 
that at sonic speed the lift-cirrve slope at zero angle of attack of air- 
foils of a single family varies inversely as the cube root of the thick- 
ness ratio. Note that as the thickness ratio goes to zero, the value of 
the lift-curve slope at zero angle of attack is indicated to become 
Infinite just as is indicated by equation (46) to be the case according 
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to lineajT theory. On the other hand, for very large a. It is plausible 
that the thickness ratio does not have any effect on c^j therefore c^ 


is proportional to the two-thirds power of a for leirge a. If hypothe- 
sis (e) is acceptable, that is, if it can be assumed that the detached 
bow wave in front of a symmetrical airfoil at an infinitesimal angle of 
attack is a normal shock wave so that (dCp/dMo) is given by equa- 


tion (27), the lift— curve slope at zero angle of attack is stationary 
with Mach number at Mq = 1 , since 


■ d 
dM^ 




Mo=i 



= 0 m 




=0 

o 


At present, the only available solution of the nonlinear equations 
for transonic flow about lifting airfoils is that of Guderley and Yoshihara 
(reference 27) for the case of sonic flow about a symmetrical double- 
wedge airfoil at an infinitesimal angle of attack. They found that 


(cj/a) = 3.32 
0=0 


(50) 


The foregoing results are summarized in graphical form in sketch (h). 
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of vanishing aspect ratio .- Two well-known results of linear 
theory are that the lift-curve slopes of wings of finite aspect ratio 
remain finite throughout the entire Mach number range and that the 
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lift-curve slopes of wings of vanishing aspect ratio are independent of 
Mach number. Therefore, eq,uation (45) implies that (CL/a)^ must he 
proportional to either for wings of vanishing A, in flow, or 
for any wing in a flow of vanishing IqC 

(C^)^ = (Cl/o.)^ = a X const. (51) 

A->o ^o"^ 


The value of the constant must be determined for each plan form by 
actually solving the equations of llnesir theory. For wings having trail- 
ing edges which possess no cutouts extending forward of the most forward 
station of ma xim uTn span, that is, triangular, rectangular, elliptical, 
etc., as well as certain swept-back wings, the value of the constant 
is n/2. (References 21, 22, 23, and 24 should be consulted for further 
discussion of this point as well as for the values of the constant for 
wings having cutouts in the trailing edge which violate the above stated 
condition.) 


It is seen from equation (5I) that the lift-curve slopes of wings 
in sonic flow decrease continuously in magnitude as the aspect ratio 
diminishes towards zero. Since, in addition, the lift-curve slope given 
by linear theory has its maximum value at Mq = 1, it is conjectured that 
the lift results of linear theory are a good approximation to those of 
nonlinear theory not only for all wings at Mach numbers far from unity 
but also for all Mach numbers for wings of sufficiently small aspect 
ratio . 


Cl 

a 


Wings of finite aspect ratio .- At the present time no solutions of 
the nonlinear theory are available for wings of finite aspect ratio. 
However, from the remarks of the pre- 
ceding pareigraphs, it is apparent that 
a curv^ representing the variation 
of C^/a with A for constant |q 
and a would have the^ollowing asymp- 
totic properties; Cl/<x would increase 
linearly with A for small A^ and be 
Independent of A for large A. In 
order to give a better idea of the 
numerical values to be expected, a set 
of tjrplcal results of this type are 
shown in sketch (l) for wings having 
triangular plan forms and symmetricaQ. 
double-wedge airfoil sections. The 
supersonic results are those of 
Stewart, Brown (references 28 and 29) >. 
and others. The subsonic results are 
those calculated by De Young and Harper 
(reference 30) using Weissinger's 
modified lifting line theory. 
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An interesting result of the foregoing remarks concerns the 
influence of thickness ratio on the lift-curve slope at zero angle of 
attack of wings in sonic flow. For wings of large aspect ratio, the 
lift-curve slope is inversely proportional to the cube root of the thick- 
ness ratio. For wings of small aspect ratio, the lift-curve slope is 
independent of the thickness ratio. 


Pitching Moment 

The remarks on the pitching-moment characteristics of 'vrLngs foUoi^ 
in a manner exactly analogous to those just stated for the lift charac- 
teristics. The corresponding statements for the pitching-moment coef- 
ficient Cm may be obtained by sin5>ly replacing Cl with Cm and L 
with M. Thus, the similaxity mles for Cm corresponding to equa- 
tions (^3) and (k 5 ) are the following, respectively: 

(C^) = [(7+1) (t/c)]^/® (Cm/a) =M(lo, A, a) (52) 



( 53 ) 


where once more M and are different functions for subsonic and 
supersonic flow. The only difference between the discussion 
of Cm and Cl is that the values of the constants of eq,uations (^), 
(50), and (51) are, of coiirse, different. Graphs of theoretical 
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pitching-moment characteristics for airfoils and for triangular wings 
corresponding to the lift results of sketches (h) and (i) are shorn in 
sketches (j) and (k) . The moment axis is teiken to he through the most 
fonrard point of the wing. 


Sometimes it is desired to 
present pitching-moment character- 
istics of wings in terms of center- 
of-pressure position rather than 
pitching -moment coefficient. Since 
the center-of -pressure position can 
he expressed in terms of Cm and Cl 
by 



( 5 ^) 


the resulting expression for the 
center-of-pressure position found 
through application of equations (^3) 
and (52) is 
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(55) 


The corresponding relation for linear theory is 


C.p.i(/uj 3 :) (56) 

Pressure Drag Due to Lift 


The similarity mle for the pressure drag of inclined wings having 
symmetrical airfoils is indicated hy equation (26\ to he the following 
if the dependence on a and t/c is written in the same manner as in 
equations ( 4 o) eind ( 4 l) ; 


Cd 


(7+lP^ 


D Uo, A, a) 


( 57 ) 
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The portion of the drag due to lift is therefore 

• aob = Ob - cdo ? a, Zy- D (Jo, A, 0) 

Da do. A, S) 


(58) 


Since Cp varies, in most cases, with the squaxe of a for at least 
a small range of a, surroimding zero (hypothesis (d)), it is advan- 
tageous to consider the dreig-rise ratio ACp/a^ rather them ACp alone, 
thus. 



1 (t/c)^/Q 

(7+1)^/® 


( 5 o>A,a) 


1 

[(r+i)(t/c)]Va 



/ ^ 


( 59 ) 


Therefore, drag-due-to-lift^data should be presented by plotting "^h^ 
variation ^d.th |q. A, and a .of a generalized drag -rise ratio ACp/a^ 
defined as follows (the prime on being dropped for simplicity); 

(ACp/a^ = [(r+l)(t/c) [ACp/a^] = (Eq> A, a) (60) 


The actual presentation of the res-ults of this three-parameter system 
may be accon^jlished as described in the section on the lift of 'vrLngs. 

Of particular interest is the simplification resulting from presenting 
only the values found at a = 0 . The simplified similarity rule is then 

(a5A^) = Da 

a=o 

For cases where linear theory applies, the following results hold; 



1 




Da, (/KT a (62) 


Wings of infinite aspect ratio.- For wings of infinite aspect ratio, 
the functional relation of linear theory for the drag due to lift, 
equation (62), reduces to 


/Aca.\ const. 


/Ac^ \ const . 

^ Uol 


( 63 ) 
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Solutions of the equations of linear theory show that the value of the 
constant is zero for subsonic flow and four for supersonic flow about any 
symmetrical airfoil. These results are valid at Mach numbers appreciably 
less than or greater than unity hut are invalid for Mach numbers near 1. 

A similarity mile for the drag due to lift of a family of affinely 
related symmetrical airfoils that is valid throughout the Mach number 
range may be obtained from eq.uation (60) by setting A = ». 




(64) 


At a Mach number of unity, eq_uation (64), for the drag due to lift, 
reduces to the following; ^ , 



(0,a), 



Mo=i 


1 

[(7+1) (t/c) 


dA 



(65) 


Equation (65), together with hypothesis (d), indicates that, at sonic 
speed, the drag-rise ratio Ac^a^ of airfoils of a single family varies 
inversely as the cube root of the thlclmess ratio. For very large values 
of qJ, the thickness ratio cannot have any effect on Ac^j therefore, 

Acd is proportional to the five -thirds power of the angle of attack. 

If hypothesis (e) is accepted, at infinitesimal angles of attack 

is stationary ’with Mach nuniber at Mq = 1> that is. 


d 


= 0, 

d 



0;=0 


ct=o 




^0’ 


Wings of vanishing aspect ratio .- Since the drag due to lift 
calculated by means of linear theory remains finite throughout the Mach 
number range, it is assumed, as in the preceding sections on lift and 
pitching moment, that linear theory is capable of describing the draig- 
due-to-lift characteristics of wings of vanishing aspect ratio at all Mach 
numbers. Therefore the following relations stemming from equation (62) 
hold: 



(67) 


Solutions of the equations of linear theory show that the value of the 
constant is rt/4 for all wings of small |Mq^ - 1 A whose trailing 
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edges possess no cutouts extending forward of ttie most forward station 
of maYiTTTinn span (i.e., triangular, rectangular, elliptical wings, etc.). 
The quoted value of the constant corresponds to the development of the 
full "leading-edge force." It is known that this force is oftentimes 
not completely realized due to a loceil separation and subseq.uent reattach- 
ment of the flow aroiuid the leading edge. If the leading— edge force is 
nonexistent, the coiresponding value for the constant is jt/2. 

Wings of finite aspect ratio. - At the present time no drag-due-to- 
lift results have been obtained from the nonlinear theory for wings of 
either finite or infinite aspect ratio. The foregoing remarks, however, 
are s uffi cient to determine that a curve representing the variation 
of with 2 for constant 5 q and a would Increase linearly 

with X for small X (unles£ the degree of attainment of the leading- ^ 
edge force also depends on A) and become Independent of % for large A. 
The resulting curve woxild presumably have the same general appearance as 
that shown in sketch (i) for 

It may sometimes be desired to present drag-due-to-lift results in 
terms of ACp/Cx,^ or ACd/oCl rather than ACp/a^. The similarity rules 
for these guantities can be quickly deduced from the foregoing results. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., March 6, 1952 
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APPENDIX A 

RELATION BETWEEN Uq AND a* STATEMENTS 
OF THE TRANSONIC -FLOW EQUATIONS 


Equations (3)j (7), (B) , (lO) , and (12) vrere presented in the text 
as being applicable to the study of transonic, as well as subsonic and 
supersonic, flow about thin wings. These equations, repeated below as 
equations (Al) through (A5), will be referred to as the Uq statement of 
the problem since the pertirrbation velocities are taken about the free- 
stream velocity Uq. The pert^arbation potential cp is defined by 

D = -UqX + (Al) 

The differential equation is . ■ , 

(1 - Mq^) Txx + Tyy + Tzz = ^ ^ ^ 9x '?xx (A2) 

Uo 

The boui'idary conditions are 
at X = - 00 

(^x)o = ("Py)o = (^z)o = 0 (A3) 

at the wing surface 

The pressure coefficient is given by 

Cp = - — (cpx) (A5) 

• o 

In the a* statement of the transonic -flow equations (equations (I3) 
through (17) in the text), it is assumed that all velocities are only 
slightly different from the critical speed of sound a*. The perturba- 
tion potential is defined by 


T’ = -a*x + $ 


and the resiilting differential equation is 

7+1 


9 


yy 


cpt = 
^ zz 


9k 9’: 


XX 


(A6) 

(A7) 


If the perturbation analysis is carried out in a completely consistent 
manner the boundary conditions are: 
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at X = -00 

('P’x)o = W- 1) - 7 ^^ ('P’y)o = ('P’z)q = 0 (A 8 ) 

at the wing surface 

and the pressure coefficient is given by 

"k f'Px - ("P'x)o ^ (AlO) 

The relation between the Uq and the a* statements can be deter- 
mined directly in the following manner. The differential equations for 
•?’ and cp will be the same if 


|2' = Z_Li |f . (1 . 

a* Bx Uq ^x ^ 


or if 


<P 


1 ’ _ cp _ a* (1 - Mq^ ) 


r. 


7+1 


(All) 


The boundary conditions for <P* corresponding to those stated for cp in 
equations (A3) and (A^l-) are 


at X = - 00 


(qp* ) = Sf rep ) - a* (1 - Mpg) ^ _ a* (1- Mo^) 

Uq ^o 7 + 1 7+1 


'■'Po " '”y>o - ' <'Pyo=§;"Pz>„ 


at the wing scrrface 


w - §j (^) -*(i) 

■7.=n ej <7— ft ' 


z=o z=o 

Finally the pressure coefficient is given by 


S’ = s 


= -^cp= _2_rHocp»4. a* (1 - 

Uo ^ Uo La* X 7+1 




(A12) 


(AI 3 ) 


(A14) 


Comparison of the above equations with those given previously for 
the completely consistent a* analysis reveals their identity. As far 
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as obtaining tbe -ralues of Cp is concerned, therefore, the a* state- 
ment of the problem may be regeirded as a transformation of the more gen- 
eral Uq statement. The resiilts so obtained are consequently valid 
throughout the Mach number range. As is evident from equation (All), 
however, the local velocities, and consequently the local Mach numbers, 
found in the a* analysis are only correct when the free-stream Mach 
number is unity. 
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■ APPENDIX B 

DEEIIVATION OF SIMILARITY RULES 


The haslc eqmtions of lineax theory and of nonlinear theory of 
transonic flow may he summarized as follows. The differential equations 
are: 


(1 + 


The boundary conditions are: 
at X = - 00 




at the wing STirface 


0 Linear (Bl) 

7 ..± 1 ^ Nonlinear (B2) 

Uo 


=0 

Sy>o V SzA 


z=o 


S(x/c) 


I 

\c h 


(B3) 




where the geometry of the wing is given hy 

(z/c) = T f(x/c, y/h) 
The pressure coefficient is given hy 

c = --^12 

P Uq ox 


(B5) 


(B6) 


If the differential equations are now transformed into a system with 
primed quantities and the proportionality or stretching fjactors are 
denoted hy s with appropriate siibscripts such that 


x’ = y’ = Syy, 

V'(l - Mo^)’ = sp Jl - 


z* 


= s^z, 9’ = S(p(p, Fq’= BuUq 

(r + i)» = r' = spT = sp(7 + 1) 


(B7) 
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equations (B1) and (B2) become 




(1-Mo^) 


I 3^9 ' 


Sq) 


Sjp 




Linear (b8) 


8^8 3 (7+1)* 89' 8^9’ 

Nonlinear (L9) 

8 8 ^ U ’ Sx’ 

T9 ° 


Similarity is assured if 9’ satisfies tbe same differential equation 
and boundary conditions as 9 • 3ji order for tbe t^ro flows to be tbe 
same, therefore, tbe following relations must bold: 


SttS 


y^3 


- 1 - 1 
— X, ■— — J., 


‘9 






A/A» 

Linear 

(BIO) 

sp2/sr 

Nonlinear 

(BU) 


where, for linear theory, A/A’ is an arbitrary constant which can be 
equated to s^^/sp if desired. 

An immediate consequence of this transformation is that the wing 
plan forms imdergo an affine transformation such that the aspect ratios 
of wings in similar flow fields ajre related, according to both linear 
and nonlinear theory, by 


S-y 1 

A* = A = A = 
sx sp , 


I 2 

1 - Mo 

(1 - Mo^’ 


(B12) 


or by 


Ji: 




^ A’ = 


Mp 2 A 


(BI 3 ) 


Since 9* is proportional to 9, the boundary conditions at x = -00 
are automatically satisfied. The boundary conditions at the wing may be 
given in either of two forms: 


59 j ^ 
5 z»/ S2 
z* =0 


. u„T. ^ 


^ 

5 z y Sz ° 5 (x/c) 
z=o 


# 0 


(Blit) 


.5^)= 

z’ =0 


11^ = 


5(x* /c’ ) ^ 5(x’/c’) 


Xi, (BI5) 

' bV 


whence, if the two wings have the same ordinate-distribution functions, 
that is, if f’(x*/cS y’/B’) = f(x/c, y/b), the ordinate -amplitude 
parameters are related as follows: 
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T» = 


^ T = 


'a / (1 -Moy 

A’ ^ 1 - Mq2 


T =< 




5s:: T. 
=r 


r (1 - 

3/2 

7 + 1 

L 1 - Mo^J 


_(7 + 1)*_ 


or as 


7(1 - Mpg) ', y 1 - Mq^ 

A’ t’ 


7(1 - 

[(7 + 1) »T' ]^/3 


7' 1 - Mo^ 

[(7 + 1)t]1/3 


Linear (BI 6 ) 

I 

T Nonlinear (BI 7 ) 

Linear (BI 8 ) 
Nonlinear (BI 9 ) 


The relation "between the pressure coefficients at corresponding 
points on the wing surface is given by 


r ' = - 


2 f Sqp»^ ^ aq) Sq)\ ^ 

Uq’ \ ^xV SyS^.V Uq 


®U®x 


Cp - < 


z‘=o 


— s 

A« ^ 


z=c 


Linear 




7+1 


(7+1) »-* 


1/3 




(B20) 
Nonlinear (B2l) 


or more completely by 


Cp' zi 

f'O 


Linear (B22) 


7+1 

LTt^ 


Nonlinear(B23) 


The foregoing relationships may be summarized in the following 
statement: The similarity rule for the pressure coefficients on a 

family of wings having their geometry given by 


(B2h) 


(z/c) =Tf(x/c, y/b) 



MCA TN 2726 


37 


is 



Linear (B25) 


^ 2/3 

n = J-L p 

^ (7+1) 1/3 


a/ 1-Mq^ 

I ( 7 + 1 ) t 1 1 / 3 ^ 


7i-Mo^A; 

C D J 


Uonlineai* (B26) 


The similarity rtiles for the lift, pitching-moment, and drag coefficients 
given hy linear theory are 


A ) linear 

*^1 ° ^ ■* ) Linear 

) Linear 

The corresponding similarity rules given by nonlinear theory are 


(B27) 

(B28) 

(B29) 


_ t ^ 3 f 5 'I 

L (7+1) 1/3 L ^ ( 7+1) t] 1/3^ V 1 "Mq AJ- 

' (^75 M 
CD=-(y^D 


nonlinear (B 30 ) 


Nonlinear (B 3 I) 


Nonlinear (B 32 ) 


It should be noted that the foregoing equations have been ‘written for 
subsoiiic flow where Mg'S 1- If 1, the radical J 1-Mq^ should 

be replaced with */Mq 3 - l. 

In the linear -theory analysis, ~K has remained a conpletely arbitrary 
coefficient to be selected as best s'ui'bs the particular problem at hand. 

For Instance, the compressible-flow relationships between two wings ha-ving 
identical pressure distributions eire found by setting A = 1. If, on -fehe 

other hand, A is set equal to */i-Mq^, the thickness ratio, camber, and 
angle of attack of related wings are identical. The greatest simplifi- 
cation of the similarity rules of linear theory occurs when A is set equal 

to /v/l-Md® /t since ‘the ninnber of parameters necessary to show the 
results of linear theory is thereby decreased by one. This degree of 



38 


MCA TW 2726 


arbitrariness in the similarity rules for linear theory is in contrast 
to the case fo?* nonlinear transonic theory in -which no undetermined 
coefficient like A is to he found. 


For the present puip)ose of gaining a better understanding of tran- 
sonic flows, the most ad-vantageous choice for A is 


A = 


12L+J^ 

T^3 


(B33) 


because then the similarity rules for linear theory assume forms identical 
to those for nonlinear transonic theory. This is important since it 
implies that solutions of linear theory and of nonlinear transonic theory 
can be expressed as functions of the same parameters and hence can both 
be plotted on a single graph in terms of one set of parameters. The t-(TO 
theories -i-rould, of course, yield "two distinct curves on such a plot. The 
curve for lineax theory would be accepted as valid for purely subsonic 
and purely supersonic flows but may or may not be valid in the transonic 
range, as discussed pre-viously. The curve for nonlinear transonic theory 
is valid not only for transonic flows, but for subsonic and supersonic 
small per-turbation flows as well. As is e-vident from the derivation of 
the basic equations, however, the results of the nonlinear transonic 
theory should be considered to be of only equal accuracy to those of 
linear theory in the definitely subsonic and supersonic regimes, despite 
the fact that -the solutions are much more difficult to obtain. 
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APPENDIX C 
SHOCK-WAVE RELATIONS 


Similarity rules were derived in appendix B througli consideration 
of the differential equation for the perturbation velocity potential cp. 
Since the transonic equation involves the existence of a velocity poten- 
tial, the derived similarity rules might be assumed limited to regions 
of flow lying between discontinuities or shock waves. It will be shown 
in the foUoiri-ng, however, that the same basic parameters govern the 
transition through weak normal or oblique shocks so that the similarity 
rules can be used to relate flows containing shock tra.ves. 

If the velocity immediately before the shock is designated by 
and the velocity conponents immediately behind the shock extending in 
directions parallel and pe ipendicular to % are designated, respec- 

respectively, by and + Ws^, the classical equation for the 

shock polar provides that 



W2^ 



U1U2 - 


— Ui^ - Uiife + 

7+1 


(Cl) 


Except for the important case of the bow wave in supersonic flow, 
is not generally alined with the direction of the x axis, but is 
inclined eL small angle. With the resolution into components parallel to 
the axes of the coordinate system and upon carrying out a small pertur- 
bation analysis analogous to that performed in the derivation of equa- 
tion (T)^ equation (Cl) provides the following relation between the 
velocity couponents (potential gradients) immediately fore and aft of the 
shock: 



This equation corresponds to the shock-polar curve for weak shock 
waves inclined at any angle between that of normal shock waves and that 
of the Mach lines. 

The striking correspondence of equations (C2) and (7) mak a it almost 
self evident that the shock relations satisfy the same similarity rules 
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as the differential equation for the perturbation velocity potential. 
This can he verified quickly hy transforming equation (C2) into a system 
1-d.th primed quantities related to the original quantities hy equa- 
tion (by), thus 



In order to assure that flows through shocks in the primed and unprimed 
systems are similar, the following relations must hold: 


Vp 


= 1, 





(C4) 


These are identical to the requirements already specified in equation (BU) 
for similarity of the flows in the shock- free regions. Therefore it 
follows that the similarity rules developed in appendix B on the basis of 
the potential equation are, in fact, valid for flows containing shock 
waves . 
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APPENDIX D 

SLOPE OF PRESSURE CURVE AT Mq = 1 


Consider a transonic flow problem in which the free -stream Mach 
number is slightly less than \mity (M5 = 1 where e is a small posi- 
tive q,uantity. Since the object is only to determine the value of a 
slope as e approaches zero, there is no loss in accuracy introduced by 
retaining only the leading powers of e. To this order, the differential 
equation for cp“ is 

The boundary conditions are: 
at X = -00 


(^x)o = ("Ppo = ^"Pz^o = ° (D2) 


at the wing surface 


(q?z) =T-. 


^oCl-S) 2=0 ^(x/c) 




!: • y) 

-hj 


The pressirre coefficient is glvep by 


^ “ ao(l-e) 


(9x) 


z=o 


(D3) 


m 


If the Mach number of the flow is now increased to Mj = 1 + e, 
keeping Sq constant, the differential equation for 9'*' is 

-26 cp^+ cpjy+ ‘Px'P^x (D5) 

hence qp+ and both satisfy the same equation provided 


'Pi + 26 cp- _26 


7+1 


(D6) 


or if 


qp+ 


1+6 

1-6 


9 


7+1 


X + const. 


(D7) 
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Consideration of tlie boundary conditions at the wing shows that the 
Tj-ing; geometry is preserved by the new potential since 


S(x/c) (l+e)aQ 


(^x) = 
z=o 


1 

(l-e)aQ 


( 9 ;) = T 
z=o 


d 

d(x/c) 


f 



(D8) 


At X = -00 , however, an alteration of the boundary conditions is 
observed since the velocity perturbations no longer vanish but take on 
a uniform value given by 


(^i)o 


4€aQ _ 4(l^-l)ao 

7+1 " ~ 7+1 ^ 


(9po = K)o = ° 


(D9) 


It can be shown that such a velocity perburbation is just that which 
would exist if there* were a nonnal shock wave of infinite lateral extent 
standing infinitely far ahead of the wing. (See eq.uation (C2).) If it 
can be assiimed that such a shock is actually present, the pressure coef- 
ficient in slightly supersonic flow is related to that in slightly sub- 
sonic flow in the following manner; 




«) = - 




2 (q>-) + 8 €_ _ ^ ^ 


ao(l-e) 7+1 


(DIO) 


z=o 


The slope of the pressure curve at Mq = 1 is therefore 


\ ^ y €— > o 

Mq = 1 


^ = 


2e 


-7+1 


(Dll) 


the value originally given by Liepmann and Bryson (reference 15) • 

It is to be en 5 )hasized that it is necessary to assume the presence 
of the velocity field at x = -00 given by equation (D9)^ due presumably 
to a normal shock wave, in both the present derivation and those of 
references 8 and 15 • V/hether or not this is a permissible assumption for 
any given case still remains an unanswered question. Intuitive consider- 
ations suggest that the results are probably applicable to symmetrical 
airfoils at^zero or infinitesimal, angles of attack but not to airfoils at 
^larger angles of attack or to wings of finite span. 
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